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ABSTRACT 

We have analyzed an XMM-Newton observation of the nearby spiral galaxy NGC 1313, which contains two 
“ultra-luminous” X-ray (ULX) sources. We measure isotropic luminosities of Lx = 2.0 x 10"^° erg u' and Lx = 
6.6 X 10^® erg u' for NGC 1313 X-1 and X-2 (0.2-10.0 keV, assuming a distance of 3.7 Mpc). The spectra 
statistically require soft and hard spectral components to describe the continuum emission; some prior studies of 
ULXs have claimed cool soft components with lower statistics. The improvement over several single-component 
models exceeds the 8cr level of confidence for X-1; the improvement for X-2 is significant at the 3cr level. The 
soft components in these ULX spectra are well-fit by multi-color disk blackbody models with color temperatures 
of kT ~ 150 eV. This temperature differs markedly from those commonly measured in the spectra of stellar-mass 
(10 Mq) black holes in their brightest states {kT ~ 1 keV). It is expected that the temperature of an accretion disk 
orbiting a black hole should decrease with increasing black hole mass. If the soft components we measure are 
due to emission from the inner region of an accretion disk, and disks extend close to the innermost stable circular 
orbit at the accretion rates being probed, the low color temperatures may be interpreted as spectroscopic evidence 
of black holes with intermediate masses: Mbh — 10^ Mq. Simple Eddington scaling arguments suggest a mimum 
mass of Mbh ^ IO^Mq. NGC 1313 X-1 and X-2 are found in optical nebulae (Pakull & Mirioni 2002), which 
may indicate that anisotropic emission geometries are unlikely to account for the fluxes observed. 


1. INTRODUCTION 

ULX sources may be defined as point-like off-nuclear X- 
ray sources in normal galaxies for which measured luminosi¬ 
ties exceed the isotropic Eddington limit for a stellar-mass 
(less than approximately 10 Mq) black hole. The existence 
of such sources was first revealed with Einstein (Eabbiano 
1989). Variability has been observed in many ULX sources 
on the timescales of months and years (in rare cases, on shorter 
timescales), which suggests that they are accreting objects. 

Although black holes with intermediate mass (IMBHs; 
10^“^ Mq) provide an attractive explanation for the nature 
of ULX sources, strong evidence for this interpretation has 
been lacking — especially from X-ray spectroscopic studies. 
Anisotropic emission from stellar-mass black holes may be able 
to account for the flux observed in some ULX sources (King 
et al. 2001). Single-component fits with the multi-color disk 
blackbody model (MCD; Mitsuda et al. 1984) have measured 
color temperatures above those commonly reported in stellar- 
mass black holes (kT = 1-2 keV; see, e.g., Sobczak et al. 2000, 
Makishima et al. 2000); as inner disk temperatures should fall 
with increasing black hole mass, these findings again point to¬ 
wards stellar-mass black holes or Kerr black holes. 

NGC 1313 is a nearby spiral galaxy {d = 3.7 Mpc, Tully 
1988). Two ULX sources — NGC 1313 X-1 X-2 — are as¬ 
sociated with this galaxy; approximately 1 kpc and 8 kpc from 
the photometric center of the galaxy, respectively (Colbert et 
al. 1995). Variability was observed in RG&477HRI observa¬ 
tions of these sources (Colbert & Ptak 2002). Colbert and 


Mushotzky (1999) claim evidence for a cool {kT =120 eV) disk 
in fits to an ASCA spectrum of NGC 1313 X-1 with a model 
consisting of MCD and power-law components; however the 
spectrum is fit acceptably by a simple power-law (T = 1.74, 
X^/d.o.f. = 244/259). Similarly, a cool disk may be implied 
in joint fits to BeppoSAXfMECS and RG&477PSPC spectra of 
M81 X-9 (La Parola et al. 2001) but not in fits to MECS and 
LECS spectra (both aboard BeppoSAX). 

We have analyzed the archival XMM-Aewfo« spectra of NGC 
1313 X-1 and X-2. Separate soft and hard components are sta¬ 
tistically required to describe the spectra. We measure low tem¬ 
peratures {kT ~ 150 eV) for the soft components. We interpret 
the soft component as arising from the inner region of an accre¬ 
tion disk, and explore the implications of cool accretion disks 
for the mass of black holes, if such objects power X-1 and X-2. 

2. DATA REDUCTION AND ANALYSIS 

NGC 1313 was observed by XMM-Newton on 17 October 
2000 starting at 03:59:23 (UT). We used only the EPIC data for 
this analysis. The EPIC cameras were operated in “PrimePull- 
Window” mode with the “medium” optical blocking filter. The 
XMM-Newton reduction and analysis suite SAS version 5.3.3 
was used to filter the standard pipeline event lists, to detect 
sources within the field, and to make spectra and responses. 
The pipeline processing and our own both failed to produce an 
event list for the pn camera. We therefore restricted our analysis 
to the MOS-1 and MOS-2 cameras. Application of the standard 
time filtering resulted in a net exposure of 29.3 ksec. 
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The source locations were determined by running the SAS 
tool “edetect_chain”. With this tool, we find NGC 1313 X- 
1 at 3fil8ml9.99s, - 66°29'10.97", and NGC 1313 X-2 at 
3hlSm2234s, -66°36'03.68" (J2000). The tool returns an er¬ 
ror of 0.2” for these positions; an uncertainty of 4” may be more 
appropriate (e.g., Foschini et al. 2002). Source counts were ex¬ 
tracted in a circle within 24” of the detected source position. 
Background counts were extracted in an annulus between 24”- 
30”. To create spectra, we then applied the selection criteria 
described in the MPE “cookbook.” These selections are as fol¬ 
lows: we set “FLAG=0” to reject events from bad pixels and 
events too close to the CCD chip edges, event patterns 0-12 
were allowed, and the MOS spectral channels were grouped by 
a factor of 15. Response files were made using the SAS tools 
“rmfgen” and “arfgen.” Spectral files were grouped to require 
at least 20 counts per bin before fitting to ensure the validity of 
statistics. 

Model spectra were fit to the data using XSPEC version 11.2 
(Arnaud 1996). The MOS-1 and MOS-2 spectra were fit jointly 
with an overall normalizing constant. The constant indicates 
that the overall flux normalizations of these cameras differ by 
less than 5%. Models were fit to the spectra in the 0.2-10.0 keV 
band. Systematic errors were not added to the spectra. Errors 
quoted in this work are at the 90% confidence level. Using the 
SAS tool “epatplot” we found photon pile-up to be negligible 
in our source spectra; this was confirmed with the HEASARC 
tool “PIMMS” using the parameters reported by Colbert & 
Mushotzky (1999) for X-1 and X-2. 

3. RESULTS 

Lightcurves of the source event lists do not show strong vari¬ 
ability on the timescale of this observation. We therefore pro¬ 
ceeded to make fits to the time-averaged spectra. The results of 
joint fits to the MOS-1 and MOS-2 spectra in the 0.2-10.0 keV 
band are listed in Table 1. 

We began by fitting single-component models commonly ap¬ 
plied to ULX sources, modified by photoelectric absorption 
(via the “phabs” model within XSPEC). For both X-1 and X- 
2, MCD, thermal Bremsstrahlung, and Raymond-Smith plasma 
models all fail to yield acceptable fits. Simple power-law mod¬ 
els provide improved but statistically unacceptable fits. Broken 
power-law models with breaks in the 0.2-10.0 keV range also 
yield improved but statistically unacceptable fits. 

An observation of NGC 1313 with Chandra finds that the 
diffuse emission in our XMM-Newton extraction region is not 
well-described with any thermal model; the best model is 
a r = 1.7 power-law with an unabsorbed flux of Tb. 3 - 7.0 = 
5.9 X 10“^'* erg/s (A. Kong, priv. comm.). This signals that 
the soft components may indeed arise in accretion disks, and 
we explored fits with a model consisting of MCD and power- 
law components, and the “CompTT” model. The MCD model 
describes a standard Shakura-Sunyaev (1973) accretion disk 
as a series of blackbody annuli. The latter model describes 
Compton-upscattering in a corona of optical depth t with elec¬ 
tron temperature kTcowna from a Wien distribution of soft seed 
photons (Titarchuk 1994); we used the version of the model 
which assumes a disk geometry for the seed photons. The coro¬ 
nal temperature could not be constrained and we therefore fixed 
kTcorona = 50 keV, which is a moderate value. X-2 is fit accept¬ 
ably by both of these models; the MCD plus power-law model 
is a significantly better fit to the spectrum of NGC X-1. The 
inability of single-component models to describe the spectra 
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Fig. 1.— The data/model ratios from fits to the MOS-1 (hlack) and 
MOS-2 (red) spectra of NGC 1313 X-1 with standard ULX models 
(the spectra are rehinned for visual clarity). The ratios correspond to 
fits with: a simple power-law, the multicolor disk hlackbody model, 
and a model with both components (see Table 1). 

of X-1 is demonstrated in Figure 1; MCD plus power-law fits 
to X-1 and X-2 are shown in Figures 2 and 3. The F-test indi¬ 
cates that the improved fits given by the MCD plus power-law 
model are significant over single-component models at more 
than the 8 ct level of confidence for X-1, and at the 3cr level 
of confidence for X-2. With the MCD plus power-law model, 
we measure isotropic luminosities of Lx-\ = 2.0;!;q3 x 10"^*^ erg/s 
and Lx -2 = 6 . 6 ;!; 2 o x 10^® erg/s (0.2-10.0 keV, assuming d = 
3.7 Mpc). If ULX sources are like Galactic black holes and 
AGN, a corona may be the source of hard X-rays. Extrapolating 
to the 0.05-100.0 keV band, we find that isotropic luminosities 
increase by approximately 50%. 

We note that some HII regions within NGC 1313 may have 
low metallicity (50-75% relative to solar abundances; Zaritsky, 
Kennicutt, & Huchra 1994). The metallicity assumed for ab¬ 
sorbing material can impact low-energy spectral measurements. 
For X-1, our best-fit cool disk plus power-law model represents 
an improvement at the 5cr level of confidence relative to simple 
single-component fits with the abundances of O and Fe fixed at 
0.25 X solar in the absorption model. With such an absorption 
model, single-component fits to the spectra of X-2 are accept¬ 
able. The low-abundance absorption model is likely inappro¬ 
priate, however; fits with this model under-estimate the depth 
of the O edge in the X-1 spectra. The Galactic column den¬ 
sity along the line of sight to NGC 1313 is 4.0 x 10^° cm“^; 
our significantly higher measurements may indicate significant 
neutral absorption within NGC 1313 or intrinsic to the sources 
(see Table 1). 

The most remarkable result of our analysis is that both mod¬ 
els require cool accretion disks {kT cf. 150 eV). These tempera¬ 
tures differ from those derived from the spectra of stellar-mass 
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Fig. 2.— The unfolded MOS-1 and MOS-2 spectra of NGC 1313 
X-1. The total spectrum, cool {kT ~ 150 eV) disk component, and 
power-law components are shown in hlack, blue, and red, respectively. 

black holes in bright states; in such spectra kT ~ 1 keV is 
typical (e.g., Sobczak et al. 2000) but temperatures can reach 
to nearly 2 keV in extreme cases. These temperatures also dif¬ 
fer considerably from those measured with single-component 
fits with the MCD model to the spectra of some ULXs, which 
approach kT ~ 2 keV (e.g., Makishima et al. 2000). 

The MCD model is based on the following relation¬ 
ship; T oc This can be used for scaling by writing 

{Mulx/Mwmq) cx {kTioM^/kTuixf. If we assume that kT = 
1 keV is a typical inner disk color temperature for 10 Mq black 
holes for values of Lx jLEdd. similar to those being observed in 
ULX sources, the best-fit MCD color temperatures measured 
for X-1 and X-2 then imply black hole primaries with masses 
near Mbh ~ 2 X 10"^ ATq. If we assume kT = 0.5 keV may be 
more appropriate for 10 Mq black holes at high accretion rates, 
this scaling still suggests masses of Mbh ~ 1.2 X 10^ Mq for 
X-1 and X-2. 

Shimura & Takahara (1995) derived a hardening factor cor¬ 
rection to account for the effects of opacity on the measured 
disk temperature and inner radius. This correction is very sim¬ 
ple; kTcorr. = r^kTobs., and ^in^corr. = vfRobs. (where / = 1.7 
is the hardening factor, and rj = 0.63 is valid for i < 70° and 
accounts for the difference between the innermost radius and 
the radius of peak temperature; see Sobczak et al. 2000 and 
Makishima et al. 2000). Merloni, Fabian, & Ross (2000) have 
found that hardening corrections may not be constant against 
changes in the physical parameters of the disk. 

The normalization of the MCD model allows for more di¬ 
rect estimates of the black hole masses if we assume that 
Rm = Risco = 8.85 km (Mbh/Mq) (where Risco is the ra¬ 
dius of the innermost stable circular orbit), appropriate for 
Schwarzschild holes. Simply, Mbh = 77/^(X/cos[i])'/^ x 
{d/lO kpc) X (8.85 km)“', where K is the model normalization, 
and i is the inclination of the system. Using the 90% confi¬ 
dence lower limit MCD normalizations (see Table 1) we find 
lower-limit masses of M^-i ^ 2200 Mq and Mx -2 ^ 830 Mq 
(assuming i = 0 and d = 3.7 Mpc; note that higher values of i 
and d would increase the minimum mass estimates, as would 
significant black hole spin). 

Defining Lndd. = 1-3 (Mbh/Mq) x 10^^ erg/s (Frank, King, 
& Raine 2002) and taking Lx in the 0.05-100.0 keV range 


Fig. 3.— The unfolded MOS-1 and MOS-2 spectra of NGC 1313 
X-2. The total spectrum, cool (kT ~ 160 eV) disk component, and 
power-law components are shown in black, blue, and red, respectively. 

as a better approximation to a bolometric luminosity, we find 
lower-limit isotropic luminosity masses of Mx-i ^ 230 Mq and 
Mx -2 ^ 70 Mq. These mass estimates are roughly an order 
of magnitude below those obtained by scaling from the MCD 
fit parameters; it is possible that X-1 and X-2 are observed at 
0.1 LEdd .7 similar to many Galactic black holes. We caution 
that the mass limits are only as good as the MCD model and 
our best-fit continuum models. 

4. DISCUSSION 

We have analyzed the EPIC MOS-1 and MOS-2 spectra of 
the ULX sources NGC 1313 X-1 and X-2. The spectra statisti¬ 
cally require soft and hard components to describe the contin¬ 
uum emission. When the soft components in X-1 and X-2 are fit 
with models for accretion disks, low disk temperatures are ob¬ 
tained (kT ~ 150 eV). Scaling these temperatures and the nor¬ 
malization of the MCD model suggests that X-1 and X-2 harbor 
black holes with Mbh — 10^ Mq, or higher. Isotropic Edding¬ 
ton luminosity scaling suggests that X-1 and X-2 may harbor 
black holes with masses on the order of Mbh — 10^ Mq. It is 
possible that we have observed X-1 and X-2 at Lx — O.lLEdd.', 
however, it is not clear which scaling method is superior. 

Colbert et al. (1995) estimate the radio power at the posi¬ 
tion of X-1 to be lO'^ W/Hz at 1.4 GHz, implying an isotropic 
radio luminosity of 10^^ erg/s. The radio to X-ray luminosity 
ratio is then approximately 5 x 10“^; this makes it unlikely that 
relativistic beaming can account for the flux of X-1, because 
beaming tends to produce flat i/Tv spectra (Fossati et al. 1998). 
Pakull and Mirioni (2002) have found that X-1 lies in the center 
of a diffuse Hq, nebula with a radius of approximately 240 pc, 
and a high [O I] A6300/Hq, ratio implying X-ray photoioniza¬ 
tion. An optical survey of massive young star clusters in nearby 
galaxies by Larsen (1999) reveals no cluster candidates within 
approximately 380 pc of our position for X-1. Similarly, Pakull 
& Mirioni (2002) find that X-2 lies at the center of an nebula 
with strong [Si II] and [O I] lines, implying that X-2 is acting 
on the local interstellar medium. Larsen (1999) reports no clus¬ 
ters within a few kpc of X-2. These findings suggest that X-1 
and X-2 may emit nearly isotropically and illuminate their local 
nebulae. It is unlikely, then, that models based on anisotropic 
emission from stellar-mass black holes (e.g.. King et al. 2001) 
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can explain the observed fluxes. 

IMBHs may be the endpoints of very massive low-metallicity 
stars (Heger et al. 2002), or perhaps Population III stars from 
the era of galaxy formation (Madau & Rees 2001). Miller and 
Hamilton (2002) have suggested that IMBHs may grow in glob¬ 
ular clusters. Ebisuzaki et al. (2001) have proposed that in¬ 
termediate mass black holes may form in young compact star 
clusters. If the disk temperature and normalization scalings are 
correct, our estimates for the masses of X-1 and X-2 are incon¬ 
sistent with black holes in X-1 and X-2 being the endpoints of 
low-metallicity or Population III stars; X-1 and X-2 may be the 
result of growth by mergers. 

Though a power-law produced a statistically acceptable fit to 
the ASCA spectrum of X-1 (x^/d.o.f. = 244/259), the values 
we have measured with the MCD plus power-law model are 
broadly consistent with those reported by Colbert & Mushotzky 
(1999) using the same model. We measure a lower disk tem¬ 
perature and harder power-law index in X-2. Makishima et al. 
(2000) fit two ASCA spectra of NGC 1313 X-2 separated by 
two years with the MCD model, and found kT = 1 .47 keV and 


kT = 1 .07 keV. Our fits to the XMM-Newton spectra of X-2 with 
only an MCD component are not acceptable, but we measure 
kT = 0.91 keV — consistent with the latter A5CA result. 

Results from a Chandra observation of NGC 5408 X-1 may 
also reveal a cool disk in a two-component spectrum (Kaaret 
et al. 2002). We speculate that future observations of ULX 
sources may also reveal cool accretion disks in some cases. Fu¬ 
ture theoretical studies of anisotropic accretion flows at high 
accretion rates may be able to explain such results in terms of 
stellar-mass black holes. 
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Spectral Fit Parameters 


Table 1 


Model/Parameter 

NGC 1313 X-1 

NGC 1313 X-2 

power-law 

Nh (102‘ cm-2) 

r 

1.90 ±0.05 

^•^0.2 

2.4±0.1 

Norm. (10^"*) 

X^Mof 

497.7/384 

T g+fl.4 

177.7/167 

MCD 

Nh (lO^* cm-2) 

0 41+0-09 

0 6+0-3 
++■'+-0.1 

kT (keV) 

1.41 ±0.06 

0 91+006 

64+P-'^^ 

Norm. (10^2) 

3 0+®'^ 

-0.4 

X^Mof 

809.5/384 

265.5/167 

Bremsstrahlung 

Nh (10^* cmT^) 

1.27 ±0.09 

1.5 ±0.02 

kT (keV) 

c q±0.6 

2.6 ±0.3 

Norm. (10-"*) 

6.4 ±0.2 

3.7 ±0.3 

X^/dof 

573.3/384 

205.6/167 

Raymond-Smith 

^H (1021 ^„-2) 

kT (keV) 

1.16±0.08 

5.4±0.3 

1.0±0.1 

T -7+0.3 

Norm. (10^^) 

1.58±0.05 

0.68 ±0.03 

X^/dof 

749.8/384 

252.2/167 

CompTT 

Nh (lO^* cm-2) 
kTseed (keV) 
kTcorona (keV) 

T 

Norm. (10^^) 

1 1+0-3 
+ ■+-0.2 
0.18±0.01 

501' 

0-80!®:°^ 

2-8Uj 

1 7+®-3 

+ •'-0.9 

0 15+®-®3 
++•+■+-0.04 

50t 

0 37+®-®’' 

++--+'-0.06 

1.7 ±0.4 

X^/dof 

431.8/383 

164.2/166 

MCD + power-law 

Nh (lO^i cm-2) 

4 4+1-2 
;->®2 

iOOU_2gA 

1 Or)±0.0o 
^•°^-0.09 

385.8/382 

12+f 

0.33 

2 0+^-^ 

^•^-0.3 
o +j±4.8 

’^•-^-0.3_ 

3+3 

-1 

kT (keV) 

Norm. 

r 

Norm. (10^) 

X^ldof 

F (10^*^ erg cnr^ 

Npower—law / Nlotal 

Lo.2-10 (10'*® erg i-')* 
io.05-100 (10"^® erg 

0-16:":*® 

1801 
-) 0+0.2 
^•■^-0.1 

3 6+'-® 

-+■ -0.8 

167.1/165 

4 0*-* 

^■++-1.2 

0.63 

0 66+®-*3 
u.oO-0 20 

1 0+®-® 

+ -++-0.I_ 


Note.— Results of fitting simple models to the EPIC MOS spectra of NGC 1313 X-1 and NGC 1313 X-2. The XSPEC model “phabs”was used to measure the 
equivalent neutral hydrogen column density along the line of sight. 

t The temperature of the up-scattering corona was not constrained by fits with the XSPEC model “CompTT” and was therefore fixed at a reasonable value. 

" The absorption-corrected or “unabsorbed” flux. 

^ The luminosity in the 0.2-10.0 keV or 0.05-100.0 keV range, assuming a distance of 3.7 Mpc. 



